Dynamical tunnelling with Bose-Einstein condensates in magnetic microtraps 
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The study of dynamical tunnelling in a periodically driven anharmonic potential probes the quantum-classical 
transition via the experimental control of the effective Planck's constant for the system. In this paper we consider 
the prospects for observing dynamical tunnelling with Bose-Einstein condensates in magnetic microtraps on 
atom chips. We outline the driven anharmonic potentials that are possible using standard magnetic traps, and 
find the Floquet spectrum for one of these as a function of the potential strength, modulation, and effective 
Planck's constant. We develop an integrable approximation to the non-integrable Hamiltonian and find that it 
can explain the behaviour of the tunnelling rate as a function of the effective Planck's constant in the regular 
region of parameter space. In the chaotic region we compare our results with the predictions of models that 
describe chaos-assisted tunnelling. Finally we examine the practicality of performing these experiments in the 
laboratory. 
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I. INTRODUCTION 

One of the most significant features of dilute gas Bose- 
Einstein condensates (BECs) is their macroscopic coherence, 
which leads to their ability to display well-known textbook 
examples of quantum mechanics on a macroscopic scale. Ex- 
amples include matter-wave interference j]]] and quantum tun- 
nelling [2]. A familiar form of quantum tunnelling occurs in 
a time-independent spatially symmetric double-well potential. 
A quantum particle initially located in one of the wells with an 
energy below the maximum of the potential barrier between 
the wells will tunnel between them, despite this being clas- 
sically forbidden (see, for instance, Ref. yfl). However, tun- 
nelling is a more general phenomenon that can be observed in 
situations involving other types of symmetries and barriers. 

In this paper we are interested in non-integrable systems, 
i.e. Hamiltonian systems with more degrees of freedom than 
constants of motion. The classical dynamics of such sys- 
tems is known to exhibit chaotic features J3]. However, for 
a Hamiltonian of the form Hq + eV where Hq is integrable 
and eV is a sufficiently small chaos-inducing perturbation (a 
quasi-integrable system), some constants of motion are locally 
conserved leading to the formation of so-called "Kolmogorov- 
Arnol'd-Moser (KAM) tori" Within these tori the dynam- 
ics of the system is still regular. Poincare surfaces of sec- 
tions of such systems show characteristic patterns of seas of 
chaos surrounding islands of regular motion [see for instance 
Fig. Ed)]. 

While particles cannot classically escape the regions bound 
by the KAM tori, quantum mechanical particles are able to 
tunnel through KAM barriers to symmetry -related islands JH]. 



Due to the similarities with spatial tunnelling in double-wells 
and the dynamical origin of the barrier the quantum particle 
crosses, this phenomenon was named dynamical tunnelling 
by Davis and Heller @]. 

A common theoretical model for dynamical tunnelling 
is that of the nonlinear pendulum, with the dimensionless 
Hamiltonian 



H= y + /<! + 2ecosf)sin 2 (|), 
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where p, q, and t are the dimensionless momentum, position, 
and time variables respectively, k is the potential strength, and 
e is the strength of the driving. This potential can be re- 
alised experimentally with cold atoms placed in an intensity 
modulated optical standing wave. In this context the obser- 
vation of dynamical tunnelling has been reported by Steck et 
al. 10] and Hensinger et al. flU]. The former used velocity- 
selected cold atoms from a magneto-optical trap, whereas the 
latter used a Bose-Einstein condensate for their initial sample 
of atoms. To reduce the effect of inter-atomic interactions in a 
BEC in the latter experiment, the BEC was allowed to expand 
until the interaction energy was negligible and as such essen- 
tially single-particle dynamics were observed. Another recent 
experiment succeeded in extracting the full phase-space rep- 
resentation of the quantum state of a kicked top undergoing 
dynamical tunnelling Jgl [Hi ■ 

One of the motivations for these studies is to explore the 
boundary between classical and quantum dynamics. The rele- 
vant parameter here is the effective Planck's constant, defined 
by the commutator [q,p] = in e s of the dimensionless posi- 
tion q and momentum p variables of the Hamiltonian, Eq. (Q]). 
By varying experimental parameters it is in principle feasible 
to carry out experiments that are identical apart from a differ- 
ent value of ^ e ff. This would, for example, lead to a variation 
in the rate of tunnelling, and this has been studied in a large 
number of theoretical papers ITll - l2lll . However, the experi- 
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merits that have been performed have only considered single, 
very large values of /z e ff, and while impressive, they have not 
directly contributed to the study of the quantum-classical tran- 
sition. 

Here we revisit the problem of observing dynamical tun- 
nelling in a experiment with ultra-cold atoms. One of the 
difficulties of working with a standing-wave potential as in 
Ref. is that it does not directly correspond to the classi- 
cal non-linear pendulum. The phase space of the pendulum 
has periodic boundary conditions for the position variable, 
while atoms in an optical standing wave can move between 
lattice sites. The value of /z e u used in Ref. iU was sufficiently 
large that after loading into the standing wave the condensate 
wave function was not split into a number of identical copies 
across multiple wells of the standing wave, but instead over- 
flowed into neighbouring wells. Thus the extent of the Husimi 
function of the tunnelling Floquet states was not confined to a 
single lattice site and long range coherence played an impor- 
tant role in the observed dynamics l22ll . This introduces some 
questions about the exact nature of the observed tunnelling as 
it does not fit the simple picture of a classical particle in a 
single well tunnelling from one period-one island to another. 
One way of avoiding this problem would be to use a trapping 
potential with a single minimum rather than multiple minima. 
Then, within the approximation that all atoms share the same 
single-particle wave function, all atoms will experience the 
same dynamics. 

For a one dimensional system, a sinusoidal modulation of 
the trapping potential can lead to a mixed phase space with 
symmetric islands of regular motion. If the potential is an- 
harmonic then dynamical tunnelling may occur between the 
period-one islands in the phase space. However, to a first ap- 
proximation almost all BEC experiments occur in harmonic 
trapping potentials. One possible realisation of an anharmonic 
potential for BEC is the radial trapping potential formed by 
the magnetic field from current-carrying microscopic wires 
combined with a homogeneous bias field on an atom chip 
rf23ll . While such potentials are harmonic at the centre, be- 
yond a certain length scale they become linear. Also, these 
microtraps can be made very tight, potentially giving access 
to a large range of h e ff , a crucial prerequisite for studies of the 
quantum-classical transition and so far not achieved with cold 
atoms. 

While quantum chaos and dynamical tunnelling are purely 
single-particle effects, the preparation of a BEC greatly sim- 
plifies the efficient loading of small regions of phase space. 
Therefore, in this paper we consider the possibilities for stud- 
ies of quantum chaos using BECs confined by magnetic po- 
tentials on an atom chip. We outline our model in section 
Hfl before describing our numerical results for dynamical tun- 
nelling in a ID atom chip potential in section [HI] We analyse 
these results using different theoretical methods, and in par- 
ticular derive an integrable approximation that explains some 
of the observed features. In section [IV] we consider the prac- 
ticalities of realising these experiments, including the effect 
of mean-field interactions in the BEC, before concluding in 
section[V] In appendix|A]we outline the derivation of the pos- 
sible potentials realizable with an atom chip, essentially con- 



sisting of an infinite wire carrying a time-dependent current 
combined with a time-dependent bias magnetic field. In ap- 
pendix[B]we derive the reduced dimensional Gross-Pitaevskii 
equation that we use for simulations in section HVl 



II. MODEL 

A. Hamiltonian 

We consider a time-dependent potential realised by a 
current-carrying wire and a homogeneous bias magnetic field. 
Reducing the problem to one dimension and making use of 
dimensionless units (see Appendix lAb. the following dimen- 
sionless Hamiltonian may be realised: 



H m = l — +k(1 + 6fCOS/)(l + q 2 ) 1/2 . 



(2) 



Here the dimensionless momentum p, position q, and time t 
variables can be converted to SI units by multiplying by the 
quantities 



P = 



B' ' 



Bo 



1 



q B' ' r £2' 



(3) 



respectively. Here B' is the magnetic field gradient at the trap 
centre, Bo is the magnetic field offset, Q is the frequency of 
modulation of the trapping potential, and m is the mass of the 
atom. The amplitude of the driving potential is given by e, and 
the strength of the potential is 



k = 



gm F fi B B' 2 
mB £l 2 



Q 2 ' 



(4) 



where g is the Lande g-factor, nif identifies the magnetic sub- 
level of the atom, ps is the Bohr magneton, and ai r is the har- 
monic trapping frequency about the minimum of the potential 
for small amplitude oscillations. A detailed derivation of the 
Hamiltonian Eq. can be found in AppendixlAl 



B. Classical dynamics 

The parameters k and e are experimentally tunable, and 
altering them allows the investigation of various dynamical 
regimes. Example Poincare sections for two values of k and e 
are shown in Fig.Q] 

In this article we consider the dynamics associated with the 
fixed points resulting from the 1:1 resonance of the unper- 
turbed motion (described by Ho = p 2 ; /2 + k y/l + q 1 ) with the 
period of the perturbation eV = e/ccos t yjl + q 2 . As the quan- 
tity k is increased from zero, two such fixed points appear at 
the origin of the phase-space at k — 1 14111 and then move 
away from each other, as plotted in Fig. [2] 

As e is increased, the KAM tori constraining the motion to 
regular behaviour are progressively destroyed, and larger re- 
gions of chaos appear about the islands centered on the period- 
one resonances. 
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Figure 1 : Poincare sections for the classical dynamics of the atom- 
chip Hamiltonian, Eq. 10 at t = 0, illustrating the influence of the 
potential strength k and modulation amplitude e on the classical dy- 
namics. We show results for the parameters: (a) k = 1.2, e = 0.1. (b) 
k = 1.2, 6 = 0.6. (c) k = 1.8, e = 0.1. (d) k = 1.8, e = 0.6. 
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Figure 2: (Color online) Distance in phase-space between the period- 
one resonances at t = as a function of potential strength k for 
three values of modulation strength e. The black curves are deter- 
mined from the Poincare section for the full Hamiltonian Eq. l[2j- 
The cyan (gray) curves are the corresponding result derived from the 
integrable approximation to the full Hamiltonian, Eq. (|9j, as intro- 
duced in Sec. lIIIBl These give an indication of the validity of the in- 
tegrable approximation. Solid lines: e = 0.1. Dashed lines: e = 0.4. 
Dash-dot lines: e = 0.8. 



erator is unitary, it has eigenvalues of the form 

/ = exp(-2mE/h eff ), (6) 



C. Quantum dynamics 

The effective Planck's constant for this system is the com- 
mutator of the dimensionless quantum position and momen- 
tum operators 

W 1 

«eff = i[p, q] = ~ ■ (5) 

Larger values of ft e ff correspond to systems that are "more 
quantum". It is possible to experimentally tune this param- 
eter while retaining fixed values of k and e. 

In a quantum system, dynamical tunnelling will take place 
between the period-one islands of regular motion J+ and T_ 
of the classical Hamiltonian because of their time-reversal 
symmetry. This can be understood using Floquet theory, a 
formalism used to describe systems with a periodic time de- 
pendence fiHUQjklSlilS- ^ Floquet operator F de- 
scribes quantum evolution for one period of the potential mod- 
ulation, To, and its eigenstates are invariant under a Tq = 2nt 
time translation. These are similar to the eigenstates of a time- 
independent Hamiltonian, which are invariant under any time 
translation. In our system, the tunnelling states are even and 
odd superpositions of states localized on the islands 1+ and 
These are expected to be eigenstates of the system's Flo- 
quet operator F,„. As noted in the introduction, li/'even) and 
l^odd) are analogous to the ground state and first excited state 
of a double-well system respectively. Because the Floquet op- 



where E is the corresponding eigenstate's quasi-energy — the 
generalization of the notion of energy to time-periodic sys- 
tems. 



D. Numerical solution of the Floquet spectrum 

To study dynamical tunnelling in this system, we must nu- 
merically determine the Floquet spectrum and identify the 
even and odd Floquet states that are localised on the period- 
one fixed points of the Poincare section. These will be de- 
pendent on the potential strength k, modulation amplitude e 
and effective Planck's constant 7z e ff ■ As the trapping potential 
is symmetric, the even and odd Floquet states are uncoupled. 
Our numerical procedure is to choose a regular grid in po- 
sition space, and then separately evolve a basis of even and 
odd combinations of position eigenfunctions for one period of 
the modulation. This determines the Floquet matrix F m for 
the Hamiltonian H m , which we diagonalize to find the Flo- 
quet states. The tunnelling states are identified as being the 
even (odd) Floquet states having the largest overlaps with an 
even (odd) superposition of coherent states centered on the 
period-one fixed points of the classical phase space. We can 
therefore expect a localised state on one island to be formed 
by \<p + ) = (|iAeven) + [lAodd))/ V2, and this will tunnel to a 
state \4>-) = (li/'even) ~ l<Aodd))/ V2 localised on the opposite 
island. The corresponding tunnelling period is given by 
^tunnel = 2nti e ff/AE, where AE is the quasi-energy splitting 
of the tunnelling states: AE = \E even - E dd\- 
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Figure 3: Quasi-energy splitting of the tunnelling Floquet states AE 
as a function of the inverse effective Planck's constant 1 /fi cff and the 
modulation strength e for the atom-chip Hamiltonian Eq. ® for two 
different potential strengths, (a) k = 1.2; (b) k = 2.0. 



Figure 4: Quasi-energy splitting of the tunnelling Floquet states AE 
as a function of the inverse effective Planck's constant l/ft e ff and the 
modulation strength e for the (a) Quartic oscillator with Hamiltonian 
given by Eq. and k = 2; (b) Nonlinear pendulum with Hamilto- 
nian given by Eq. l[8) and k = 1.5. 



III. RESULTS 

In this section we investigate the dependence of the quasi- 
energy splitting AE on the tunable parameters for our system 
k, e, and tt e g. Typical results are shown in Fig.[3]for (a) k — 1 .2 
and (b) k = 2.0. These figures have several noteworthy fea- 
tures, many of which have been observed previously in other 
systems 03 EES Si- 



A. Overview 

For our analysis, we divide the parameter space shown in 
Fig.[3]into two regions according to the characteristic behav- 
ior of the quasi-energy splitting. Firstly, we see in the bottom- 
left corner of Fig. Oa-b) that for a large parameter range the 
behaviour of AE is not completely unruly, but instead some 



grooves and plateau structures are obvious. The dependence 
of AE on l//z e ff differs from the exponential behaviour char- 
acteristic of Hamiltonians of the form H — p 1 12 + V(q), while 
still being quite smooth. In this regime, which we call quan- 
tum/regular regime (QRR), the quasi-energy splitting period- 
ically falls to zero. To the best of our knowledge, this has 
not been commented on in detail in previous studies — we 
account for the origin of these grooved structures in section 
IIIIBI This result constitutes the most significnate theoretical 
development of this paper. 

Secondly, in the top right corners of Fig. [3ja-b), classi- 
cal chaos has the strongest influence on the system. This 
regime is characterised by the dramatic and apparently dis- 
orderly fluctuations of AE lfl2l [l3l OH], and we refer to it 
as classical/chaotic regime (CCR). In our discussion of the 
CRR, section HIlCI we apply existing theories by Podolskiy 
and Narimanov [19;] and Eltschka and Schlagheck 111611 to our 
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specific scenario. Comparison of these theories with our nu- 
merical simulations shows good agreement and completes our 
analysis of the present atom chip system in the single-particle 
regime. 

We note that the division described above is not peculiar to 
the atom-chip Hamiltonian, but appears to be a general feature 
of driven quasi-integrable systems. We have also diagonalized 
the Floquet operator for other nonlinear potentials such as the 
quartic oscillator 

rP- 

H = — +k(1 +ecosf)<? 4 , (7) 
as well as the nonlinear pendulum 

P 2 

H = + ecos/)cos<7, (8) 

and found the same features. Sample results are shown in 
Fig. |U The latter Hamiltonian has been studied elsewhere 
both theoretically d El HIE! and experimentally HH. 

B. Quantum/regular regime 

We define the quantum/regular regime as the region in 
parameter space where the quasi-energy splitting AE varies 



We find that this integrable approximation represents the clas- 
sical dynamics of H„, reasonably well up to k = 1.3 and 
e = 0.6. This can be seen by comparing the Poincare sections 
for the full Hamiltonian Eq. and the integrable approxi- 
mation Eq. (|9]i in Fig. [6] Also, the prediction of the distance 
between the period-one fixed points of the integrable approx- 
imation Eq. Q and the full Hamiltonian Eq. (|2]i is shown in 
Fig. |2] as a function of k for three values of e. 

To proceed we quantize the integrable Hamiltonian Hi, 
Eq. (|9). In principle we could choose any operator ordering; 
given the symmetry of Eq. (0 we choose symmetric ordering. 
This is then easily diagonalised numerically, and the resulting 
spectrum is plotted in Fig. [5jb). This should be compared to 
the spectrum of the Floquet operator F„, for the full Hamil- 
tonian shown in Fig. 13a). Qualitatively the spectra are simi- 
lar and quantitatively the prediction for the quasi-energy split- 
tings of the tunnelling states are in good agreement with that 
found from the full Hamiltonian as shown in Fig. EJc). We 
conclude that the tunnelling states of the integrable approx- 
imation contain the essential features of the full system for 
k < 1.3, e < 0.6, and 1 > 1. It was previously shown that 
the enhancement of tunnelling does not require chaos but can 
originate from avoided crossings in the Floquet spectrum l27ll . 



smoothly, apart from the apparently periodic troughs where 
the splitting goes to zero. In both Fig. [3ja) and (b) it is the 
region for which e < 0.3. The origin of this behaviour can 
be seen by plotting the quasi-energy spectrum of the Floquet 
operator F m as in Fig.|5ja-b). We find that at regularly-spaced 
values of 1 /ties, a Floquet state not involved in the tunnelling 
deflects the tunnelling state of the same parity through an 
avoided crossing scheme very similar to the result of a typ- 
ical first-order perturbation theory calculation. 



This scenario has previously been invoked in the context of 
tunnelling or dynamical tunnelling suppression i25l[26ll. often 
in order to describe the phenomena inducing the CCR fluc- 
tuations Il2l El Ultl ■ More generally, it has been noted that 
level repulsion occurs for dynamical tunnelling in the pres- 
ence of a chaotic sea Oil EH EH- However, we demonstrate 
here that this effect has no chaotic origin. This conclusion was 
also reached in i27ll by different means. Here we make use of 
secular perturbation theory 10], and assume that k — 1 is of 
order e and therefore within the islands of regular motion we 
have both q — 0{ -\[e) and p = 0( -\fe). This allows us to to 
derive a second-order integrable approximation Hi to the full 
atom-chip Hamiltonian H m given by Eq. (0 as |0] 



(9) 

I 

Here we have demonstrated that even non-integrability is not 
necessary. 

To understand the quasi-periodic vanishing of the difference 
in quasi-energy of the tunnelling Floquet states as a function 
of h e s, we first plot in Fig |7Ja) the energy of the classical 
trajectories of Hi not located on the islands for k — 1.2 and 
e = 0.1 in Fig. |7ja) as a function of the inverse action J~ l , 
where 10] 

J=^(j)pdg. (10) 

A gap exists in the curve for the range of the inverse action 
2.5 < J~ l <,4.2 corresponding to energies E > 1.29 x 10~ 2 
for which the classical trajectories are instead located on the 
islands. 

We can now apply the Einstein-Brillouin-Keller (EBK) 
quantization method, which states that in order to correspond 
to a quantum eigenstate, the trajectories must satisfy 

/ = H eB (n + 1/2) with n e N. (11) 

This will discretise the curve of energy versus inverse action 
in Fig. |7{a) for any particular value of /z e ff . Now as 1 /ti e g is in- 
creased, the spacing of the points along the x-axis will also be 



H t = 



3 /, 2 \3 3(jc-2)/ 2 2 n2 6(k- 1)-3(k- l) 2 + e 2 2 6{k - 1) - 3(* - l) 2 - 5e 2 , 
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Figure 5: (Color online) Comparison of the Floquet spectrum of the full Hamiltonian H m , Eq. 10, and the energy spectrum of the integrable 
approximation Hj, Eq. (|9), for k = 1.2 and e = 0.1. (a) The quasi-energy spectrum of the full Hamiltonian as a function of l//? e ff- The black 
points are for the even Floquet states, and the cyan (gray) points are for the odd Floquet states. The points corresponding to the tunnelling 
states are joined with thick solid lines of the same color, (b) The energy spectrum of the integrable Hamiltonian as a function of 1/%b, with 
the same color scheme as for (a). The structures here are very similar to those in (a), indicating that the integrable approximation captures the 
important features of the full Hamiltonian in this regime, (c) A plot of the quasi-energy difference of the even and odd tunnelling states as a 
function of 1//W- Black curve: full Hamiltonian. Cyan (gray) curve: integrable approximation. 



increased, and particular eigenstates will move from the left- 
hand part of the curve (outside of the islands) to the right-hand 
part (inside the islands), and have a strong effect on the tun- 
nelling states. The energy of EBK quantized trajectories for 
n — to n — 10 is plotted in Fig. 13b) — even trajectories as 
thin black solid curves and odd trajectories as thin cyan (gray) 
dashed curves. The thick solid black curve is the energy of the 
even tunnelling state, and the thick solid cyan (gray) curve is 
the odd tunnelling state, both found from numerical diagonali- 
sation of the integrable approximation to the full Hamiltonian, 
Eq. (0. In Fig.[5]it is easy to see the avoided level crossings 
between Floquet states with the same parity. The curve for 
the quasi-energy of each tunnelling Floquet state changes di- 
rection at each avoided crossing, resulting in regular crossings 
of the quasi-energies of the odd and even Floquet tunnelling 
states, resulting in AE = and an infinite tunnelling period. 
It seems quite reasonable that this explanation applies directly 
to the quasi-integrable Hamiltonian (0 in the QRR regime. 

We have found that the spacing of the occurrences of AE = 
as a function of 1 //z e ff decreases with increasing k. For ex- 
ample, it can be seen in Fig.[3]that the spacing of the valleys 
as a function of l/ti^s decreases between (a) with k — 1.2 
and (b) with k = 2.0. This can be explained in the follow- 
ing manner. We find that that the period-one islands remain 
of a similar phase space area for a given k. Let the minimum 
action of any closed trajectory that lies entirely outside the is- 
lands be denoted as J\ . Likewise, let the maximum action of 
any closed trajectory that lies entirely inside the islands be de- 
noted as J2- Then the 77th (with n large enough) crossing will 
occur between l/%s = (n + 3/2)/ J\ and (n + 1/2)/ Ji- As k 
increases, the distance between the islands in phase space in- 
creases, meaning J\ and J2 increase, and the spacing between 
the nth and (« + 1 )th quasi-energy degeneracy of the tunnelling 
states will decrease. This is in agreement with the behavior of 



the Floquet spectrum of the full system. 

In summary, in this section we have shown that the quan- 
tised integrable approximation Eq. (0 accounts for all the sig- 
nificant features of dynamical tunnelling of the full system 
Hamiltonian Eq. (|2]) in the QRR. Therefore we have demon- 
strated that the underlying classical non-integrability of the 
full system has little effect in the QRR. 



C. Classical/chaotic regime 

The defining characteristic of the classical/chaotic regime is 
the irregular and large fluctuations of the quasi-energy split- 
ting AE as the parameters 1 /fi e s and e are varied. Before con- 
sidering the fluctuations, we comment on the coarse behaviour 
of the tunnel splitting as we enter the CCR from the QRR, by 
choosing e near the upper boundaries of Fig. [3ja-b) and then 
increasing 1 /^eft- 
Averaging over the fluctuations of the quasi-energy differ- 
ence of the Floquet states as a function of 1 /ti e s, we find that 
the overall tendency of the tunnelling splitting variation re- 
sembles an exponential decrease with increasing l£h e g, as ex- 
pected for dynamical tunnelling phenomena lfl7l 12811 . It is not 
surprising that the tunnel splitting tends to decrease as the po- 
tential strength k is increased and the two islands move away 
from each other in phase-space. However, it should be noted 
that after initially decreasing with increasing e, AE starts in- 
creasing again as the system enters the CCR lfl2tl . This lat- 
ter phenomenon can be attributed to so-called "chaos-assisted 
tunnelling". This occurs when the dominant transport mecha- 
nism from one island to the other is no longer direct quantum 
transport as in a regular system, but instead a three-step pro- 
cess. Firstly particles tunnel from the centre of one island to 
the sea of chaos, then are classically transported through the 
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Figure 6: Comparison of the Poincare sections for the full Hamil- 
tonian H„„ Eq. {2), and the integrable approximation ff;, Eq. ^ at 
f = 0, for (a) k = 1.3 and e = 0.1, (b) k = 1.1 and e = 0.6, (c) k = 1.3 
and e = 0.6. Although the presence of chaos for H„, in this latter 
regime makes the Poincare sections of H m and Hj look rather differ- 
ent, we find good agreement for the regular regions. The quantum 
behaviour of the two systems agree for relatively large values of 
as shown in section llllBI 



chaotic sea from the vicinity of this island to the vicinity of 
the other, and then eventually tunnel from the sea of chaos 
to the center of this latter island lfl9tl . The complete behav- 
ior of AE as 1/^eff is increased is shown in Figs. [8] and [9] A 
number of regions with distinct changes in the average slope 
of the log 10 (AE) versus l/ti e s curve are evident. These can 
indeed be partially accounted for by changes of the dominant 
transport mechanism from one island to the other. 

For large values of ft e ff, we expect the quantum behaviour 
of chaotic systems to approach that of regular systems, as 
the quantum "coarse-graining" makes them insensitive to the 
presence of fine phase-space structures. We are hence in the 
QRR, where the quasi-energy splittings of the time-modulated 
system agree with those of the integrable approximation, 
showing that the dominant transport mechanism from one is- 
land to the other is direct tunnelling. 

Starting from the QRR and decreasing the effective 
Planck's constant, we reach a regime where Ti e s is comparable 
to the phase space area of the period-one islands of regular 
motion. If chaos is present in the classical phase space, the 
dominant tunnelling mechanism is then expected to be chaos- 
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Figure 7: (Color online) (a) Energy versus inverse action for the tra- 
jectories of the integrable approximation Hj that are not located on 
the islands (solid black curve). The gap occurs due to the presence 
of the islands. The inset shows selected trajectories of Hj. The en- 
ergy of the trajectories in cyan (gray) are denoted in the main figure 
by cyan (gray) points as indicated by the arrows, (b) Energy of the 
semi-classical and quantum states as a function of 1 /S eff . The thin 
curves are the EBK quantization of the trajectories of Hj, Eq. J9). 
The thin solid black curves correspond to even parity states, and the 
thin dashed cyan (gray) curves to odd parity states. The thick black 
curve is the energy of the even tunnelling state, and the thick cyan 
(gray) curve is the energy of the odd tunnelling state, both found via 
the numerical diagonalization of Hj, Eq. (|9j. The avoided crossings 
of the quantum tunnelling states with other states of the same parity 
are clearly visible. 



assisted tunnelling as described earlier. According to Podol- 
skiy and Narimanov fl9ll . the tunnel splittings yielded by such 
a mechanism have the following dependence on h e ff: 



AE oc ti e ff 



Y(A/7Tn eS ,2A/7Tn eB ) 
r(A/^eff + i,o) ' 



(12) 



where A is the phase-space area of one island of regular mo- 
tion and F the upper incomplete Gamma function 12911 . We 
observe the transition from the QRR to this regime in Fig. [8] 
where the tunnelling behaviour of H,„ diverges from the ap- 
proximate exponential behaviour yielded by H, to follow the 
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Figure 8: (Color online) Quasi-energy splitting of the tunnelling Flo- 
quet states versus l/fi e ff for k = 1.2, e = 0.75. Thin solid black 
line: exact result. Thick cyan (gray) solid line: numerical diago- 
nalization of the integrable approximation. Dashed red line: result 
of the Podolskiy-Narimanov theory with the proportionality coeffi- 
cient of Eq. i 1 2t as a fitting parameter. Inset: Poincare sections of 
the full Hamiltonian (right) and of its integrable approximation (left). 
The tunnelling behaviour of the full Hamiltonian H m [Eq. 10] agrees 
with the integrable approximation H, [Eq. l[9j] as long as the quan- 
tum coarse-graining effect prevents the quantum particle from seeing 
the layers of chaos. For smaller values of ft cff , tunnelling is enhanced 
by the presence of chaos in the underlying classical phase-space ac- 
cording to the Podolskiy-Narimanov model fl9Tl . 




5 10 15 20 

Figure 9: (Color online) Quasi-energy splitting of the tunnelling Flo- 
quet states as a function of 7^ for k = 2, e = 15/32. Solid black line: 
exact numerical result. Red dashed line: Podolskiy-Narimanov the- 
ory (one fitting parameter). Blue dash-dot line: Eltschka-Schlagheck 
theory (no fitting parameter). Inset: Poincare section centered on 
I + . One secondary resonances chain is clearly dominant. For large 
values of the quantum coarse-graining prevents the system from 
seeing these small phase-space structures, and its behaviour is then 
well described by Eq. d!2l >. For smaller values of S cff the Eltschka- 
Schlagheck model llfll . which takes the island chain into account, 
shows much better agreement with the numerical data. 



Podolskiy-Narimanov theory as ti e $ is decreased. A more 
generic result than Eq. ( fTSl i for the chaos-assisted tunnelling 
splitting has recently been found by Backer et al. i2lll . 

However, neither the Podolskiy-Narimanov model nor the 
Backer et al. model takes into account fine features of the clas- 
sical phase space, for example partial dynamical barriers em- 
bedded in the sea of chaos DOIl . or the internal structure of 
the islands of regular motion themselves. Therefore, even if 
some systems are well described by these models for a certain 
range of /z e ff , decreasing this parameter and therefore the res- 
olution of the quantum coarse-graining can make the system 
"see" those small structures. As can be seen in Fig. [9] this 
leads to dramatic deviations from Eq. (112) . 

For ft e ff < 0.1, it is once again possible to explain the ob- 
served deviations by a change in the dominant transport mech- 
anism from one island to the other. It is known that in the 
semiclassical limit h e g — » 0, tunnelling through regular phase- 
space is greatly facilitated by the presence of resonance chains 
iflii [I5I1 . It is therefore expected that for low h e g tunnelling 
between an island of regular motion and the sea of chaos will 
be enhanced by the presence of a secondary resonance chain 
embedded in this island. 

A model reflecting this idea has been developed by Eltschka 
and Schlagheck llal . with a detailed description given in 
Ref. If20ll . Let us assume each of the symmetric islands of 



regular motion supports a r.s secondary resonance chain — 
i.e. where s internal oscillation periods match r driving pe- 
riods, and r sub-islands are visible on the Poincare section. 
This situation is well illustrated in our system for the param- 
eters presented in Fig. [9] (see inset). Firstly, the authors ap- 
proximate the dynamics within the island to that of a pen- 
dulum and treat the pendulum potential 2V ns cos rO as a per- 
turbation. This perturbation couples the ground state |i/<o) of 
the unperturbed approximate Hamiltonian only to the excited 
states \4ii r ) where / is an integer. In the following, we will de- 
note = \iffi r ). Denoting by 2nl c the phase-space area of 
the island of regular motion (measured on a Poincare section), 
they calculate the coupling of the unperturbed ground state 
of the island to its Mi 0-type state through this resonance, 
where k is defined by I{k-\) < I c < h I e N denotes the 
action associated with the /th 0-state). In other words, they 
calculate the coupling of the approximate ground state of the 
island to the lowest 0-state localised in a region where chaos 
is present for the exact Hamiltonian. They then assume that 
this coupling describes the effective coupling of <pu to the sea 
of chaos and model the latter by two matrices of the Gaussian 
orthogonal ensemble (one for each parity), which means that 
they neglect the effects of partial barriers in the chaotic part of 
phase space. Using the results from Ref. lfl8tl . they are there- 
fore able to calculate an expectation value for the logarithm of 
the quasi-energy splittings for each set of parameters, which 
we compare to our numerical results in Fig. [9] 



We find reasonable qualitative and quantitative agreement 
in the range h e g < 0.1, which is expected as this is precisely 
the order of magnitude of the phase-space area of the dom- 
inant chain's resonances. We also observe an effect already 
noted by Eltschka and Schlagheck [16], namely the fact that 
the sharp decrease of the analytically calculated tunnelling 
splittings is shifted to lower values of %s as compared to the 
similar decrease of the splittings obtained from the numerical 
simulations. Such steps occur when the value of the index k 
changes, i.e. every time an unperturbed </>-state of the approx- 
imate Hamiltonian crosses the somewhat artificial limit I c of 
the island of regular motion. This happens for every integer / 
such that I c = H e g(lr + 1 /2) (in our case, r — 6 and / = 1). It 
is known, however, that states that are neither chaotic nor reg- 
ular (hierarchical states: see Ref. i3ltl ) can be localised in the 
vicinity of the island. Therefore, the "effective value" of the 
islands' phase-space area could be a little higher than 2n x I c . 
This could account for the observed shift of the step. 

Recently a unified description that can predict tunnel split- 
tings from the direct tunnelling regime to the resonance as- 
sisted regime has been developed 1331 . It is hopeful that this 
theory would provide an improved description of the the data 
in Fig. [9] however its complexity puts it beyond the scope of 
this paper. 

IV. PROSPECTS FOR EXPERIMENTS 

We now examine the feasibility of performing dynamical 
tunnelling experiments that explore the effects presented in 
the previous sections with a Bose-Einstein condensate in a 
magnetic microtrap. To this end we discuss the realisation of 
a one-dimensional Hamiltonian in section IIV Al realistic pa- 
rameters for the trapping potential in section HV Bl the effects 
of atomic interactions in section HV CI the initial state creation 
and loading in section HV Dl and finally condensate sizes and 
inelastic losses in section HV El 



A. Confinement geometry 

A key difficulty is that the form of the Hamiltonian Eq. © 
is realized in terms of the two-dimensional (2D) radial coor- 
dinate r = -y/x 2 + y 2 of the atom-chip trap, see appendix [A] 
Instead, we require it to describe a ID degree of freedom, 
say q — x. Without this constraint, the islands of classi- 
cally bounded regular motion can be dynamically linked via 
the second radial dimension (y). Classical diffusion can then 
populate the other island and mask the signature of dynami- 
cal tunnelling. While the tunnelling and diffusion rates will 
in general be quite different, it would be desirable to avoid 
this effect altogether by freezing out the dynamics in the y- 
direction. This could be achieved by applying an optical lat- 
tice formed by two laser beams propagating along the ±y di- 
rections and confining the atoms in one well of the lattice, as 
sketched in Fig.fTOl 

In contrast to the y-dimension, the z-dimension can be very 
weakly confined on an atom-chip. The condensate in that di- 
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Figure 10: (Color online) Schematic of the arrangement of conden- 
sate, confining potentials and atom-chip surface (yellow) with char- 
acteristic z-shaped wire. The atom chip potential has the form Eq. (TJJ 
along x = q (red) and is weakly harmonic along z (potential not 
shown). Superimposed is an optical lattice potential (blue), freezing 
out the dynamics along y. The gray shaded volume is the result- 
ing isodensity surface of a trapped BEC. Embedded density profiles 
along the x, y, z, co-ordinate directions are also sketched. 



rection then has an approximately "infinite" extension. We 
will check this assumption below in section HV El 



B. Accessible tunnelling parameters 

We begin by proposing a regime in which it is feasible to 
observe dynamical tunnelling. From a practical point of view, 
we wish to minimise the tunnelling period and the sensitiv- 
ity of the tunnelling to any small fluctuations in experimental 
parameters. Distinguishable period-one islands appear in the 
Poincare section for k > 1, and the closer the islands are in 
phase space, the shorter the tunnelling period. However, if 
the two islands are too close to each other, separating them in 
time-of-flight absorption images could be challenging. High 
values of fi^ yield faster tunnelling, but require very tight 
trapping potentials. A convenient and accessible choice of 
parameters within the QRR is k — 1.01, ft e ff = 1/73, and 
e = 0.48, which gives a relatively short tunnelling period of 
159 modulations of the trapping potential. This has the ad- 
ditional advantage that the tunnelling period is relatively in- 
sensitive to small, uncontrolled experimental variations in the 
exact value of the parameters, with the strongest dependence 
being on k. We calculate that the following variations would 
induce shifts of less than 10% in the tunnelling period: 

1.009 < k < 1.011, 1/77 < 7j eff < 1/69, 0.37 < e < 0.59. 

(13) 

The Poincare section and Husimi functions of the two tun- 
nelling Floquet states are shown in Fig. QT| To achieve this 
with a BEC of 87 Rb atoms in the F = 2, nif = 2 magnetic sub- 
state confined on an atom chip would require a current through 
the trapping wire of / = 63 mA and a bias magnetic field of 
Bt — 5 G. In the limit of an infinitesimal wire, this yields a 
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Figure 11: (a) Poincare section at t = for the proposed experimen- 
tal parameters k = 1.01, fe e ff = 1/73, e = 0.48. Note that no chaos 
is present between the islands, (b) Husimi function for the even tun- 
nelling state, (c) Husimi function for the odd tunnelling state. The 
axes for (b) and (c) are the same as for the main graph (a). 



trap 25 //m from the wire which, when combined with an axial 
magnetic field of Bq = 0.25 G, has a radial trapping frequency 
in the harmonic approximation of co r = 2tt x 5 kHz. At such 
distances wire imperfections that could distort the confining 
potential and cause cloud fragmentation can be avoided lO . 
The trap depth is of the same order of magnitude as the atomic 
Zeeman energy from the magnetic field due to the wire current 
at the location of the trap minimum. For our parameters it is 
roughly equal to E m . dg = gmffisBi, ~ 0.3 mK, more than suffi- 
cient to trap a BEC at nanokelvin temperatures. 



C. Bose-Einstein condensates and mean-field effects 

So far we have considered the physics of single atoms 
within the driven potential. This has previously been sufficient 
to describe experiments in very dilute thermal or condensed 
gases JtI Htl ■ However, to work with a trapping potential with 
a single minimum, to controllably load small regions of phase 
space, and to still be able to image the resulting dynamics, it is 
necessary to use atoms that are sufficiently cold and dense that 
they are Bose condensed. In realistic BECs nonlinear mean- 
field interactions arising from atomic s-wave collisions can 
play an important role in the dynamics l34ll . These can be de- 
scribed using the one-dimensional Gross-Pitaevskii equation 
(GPE) 



„ d 
in e s—if/ = 
ot 



.-£.— + V(q,t) + U iD W 
2 oq- 



ifr. (14) 



Here the condensate wavefunction if/(q, t) is normalized to one 
and the potential is V(q,t) = k(1 + ecos(/))(l + q 2 ) 1/2 . The 
relation between the effective one-dimensional interaction- 
strength Uiu and the physical parameters of the system is de- 
rived in appendixlBl 
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Figure 12: (Color online) (a-c) Dynamical tunnelling signatures from 
condensate mean-field theory according to Eq. l !14b . for parameters 
as in Fig.[TT]and U m = lx 10"*. (a) Husimi function (FWHM) of the 
initial state located within the Poincare section: (red) initial Floquet 
state |</> + >. (blue) experimentally accessible Gaussian approximation 
as explained in the text, |i/> + ). (b) Dynamical tunnelling in momen- 
tum space from initial Floquet state and (c) from initial Gaussian 
approximation as explained in the text, (d-f) The same as (a-c) but 
for parameters k = 1.3, %k = 0.5, e = 0.2 and t/io = 0.01, with 
p = 0.65 and /q n i = 0.4. 



We have performed numerical simulations of the ID Gross- 
Pitaevskii equation ([Pit with the BEC intially loaded in the 
pure Floquet superposition state \<p+). For the case described 
in section UVBI we find that dynamical tunnelling oscillations 
can be observed with Eq. (TBI) for nonlinearities up to about 
U m - dx — 1.8 x 10~ 4 , beyond which the nonlinearity shuts down 
dynamical tunnelling lf35ll . Here we choose Um — 1 x 10 4 
resulting in unhindered dynamical tunnelling oscillations as 
demonstrated in Fig. [T2l b). where we show the momentum- 
space distribution at integer multiples of the modulation pe- 
riod. This momentum-space image constitutes a direct exper- 
imental observable, as standard time-of-fiight expansion and 
absorption imaging techniques convert the centre of mass mo- 
mentum into spatial positions that are easily resolved l36ll . An 
improved signal to noise can potentially be obtained by imag- 
ing along the weakly-trapped z-direction, as the atomic den- 
sity will be spread across fewer pixels on the camera. 



D. Experimental and numerical initial states 

In this section we consider how to best load the initial Flo- 
quet superposition experimentally, as well as the effects 
of imperfect loading. In Fig. fT2l c) we show the results of 
a simulation with the same parameters as for Fig. [T2Tb). but 
with the initial state being a BEC dynamically prepared from 
the ground state of the trapping potential such that it obtains 
a large overlap with the Floquet superposition \(f>+). We la- 
bel the resulting state \<p+). Numerically, we create \tp + ) as 
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follows: (i) We assume an initial potential strength K\ ni , differ- 
ent from k. In the potential V(q) = Aj n i(l + + q 2 ) 1 ^ 2 , we 
determine the condensate ground state using imaginary time 
evolution ll36ll . (ii) The resulting wave function f\{q) is given 
a momentum kick fi = f\ exp [ipotieffi] . (iii) We vary and 
po to maximize the overlap (0+1/2) and finally take this as the 
initial state po] for our simulation. We see the dy- 

namical tunnelling arising from this procedure in Fig. [T2l c). 
where we have used po = 0.244 and = 4.5. 

The corresponding experimental sequence is as follows: we 
begin with a stationary BEC trapped in the ground state of 
the magnetic potential with trap strength given by Kj n i. The 
wire current would then be slowly decreased by an amount 
51 = 0.78 mA and then suddenly switched back to 63 mA to 
start the cloud oscillating about the centre of the trap mini- 
mum, acquiring the maximal momentum po. When the BEC 
passes the bottom of the trap for the first time (t = 0), the 
simultaneous modulation of wire current and bias magnetic 
fields at e = 0.48 and O = 2n x 4.975 kHz begins, resulting 
in the loading of an equal superposition of the odd and even 
tunnelling Floquet states with a fidelity of up to 95%. At this 
moment the trap-strength is also changed to k. Following the 
described sequence, the state \<p+) can be directly created in 
the experiment. 

On closer inspection, the simulations beginning with this 
initial state show a slightly different dynamical tunnelling pe- 
riod compared to those that begin from a Floquet state su- 
perposition. We find such modifications in the presence of 
nonlinear interactions whenever the initial state slightly dif- 
fers from the exact Floquet state \<p+). As we vary Aj n i and po 
used to define \<p+), the tunnelling period undergoes small but 
continuous changes. This effect is absent in the linear case 
with U = 0. 



E. Condensate parameters 

The parameters of Fig. fT2l a-c) amount to a condensate of 
N - 171 atoms [Eq. dB7H , assuming the axial extension of 
the condensate is 50 pm, with a corresponding peak density 
of 2 x 10 14 crrT 3 . For comparison, we show in Fig. fT2f d-f) a 
second parameter set with much larger h e ff = 1 /2, with other 
parameters k = 1.3, e = 0.2 and U\d = 0.01, corresponding 
to N = 20 atoms. Realising 7t e ff = 1/2 would either require us 
to lower Bq to 17 mG or increase o) x to about 2k x 300 kHz. 
While challenging, precise imaging of small clouds of atoms 
as discussed here has been achieved 1371 13811 . however the re- 
liability of mean-field theory in this regime is questionable. 

We estimate that three-body losses in this scenario would 
limit the BEC lifetime to a few hundred milliseconds lf39tl 
which is sufficiently long to conduct the experiment. The 
Majorana transition-induced losses will be negligible as 
WLarmor/^7 ^ 10 at all times, where 6 is the angle between 
the magnetic field and the z axis at the position of a classical 
period-one resonance. 

As mentioned in section IIV Al the freezing of the y- 
direction could be achieved by applying an optical lattice 
formed by two laser beams propagating along the +y direc- 



tions and confining the atoms in one well of the lattice. The 
transverse (radial) size of the condensate with 171 atoms is 
less than 200 nm, which would easily fit into a single lat- 
tice site. Therefore, not only would dynamical tunnelling take 
place exclusively in the plane perpendicular to the lattice, it 
would also occur in a single lattice site. 

Finally, we have assessed the influence of the trapping po- 
tential in the z-dimension, and of experimental imperfections 
in the trap alignment. We have explicitly verified that the 
weakly confined z-dimension can be ignored, as mentioned 
in section IIV Al To this end, we modelled the dynamics of 
Fig- Housing the full two-dimensional GPE, Eq. dB3l ), with a 
weak harmonic trap in the z-direction. The axial condensate 
wave function had a non-trivial Thomas-Fermi profile l36Tl . 
The tunnelling signal in momentum-space is as clear as those 
shown in Fig. fTZt b.c). even accounting for imaging along the 
z-axis. Furthermore, in two dimensional simulations we have 
found the dynamical tunnelling to be robust to changes in the 
initial state away from the ideal Floquet superposition, as well 
as small offsets from the trap centre. 

Unfortunately, it seems that the tight trapping potentials 
combined with the small atom numbers present a significant 
experimental challenge for dynamical tunnelling with BECs 
to be realised as described here. The small atom numbers 
found for these examples are mostly due to the need to freeze 
out dynamics in the y-direction. An investigation to what ex- 
tent dynamical tunnelling can still be identified if dynamics in 
the y-direction is retained would be of interest, but is beyond 
the scope of the present paper. Note that the small number of 
atoms is required in order to keep the peak-density low enough 
for the effect of interactions to remain small. The effect of 
stronger interactions can be quite complicated, but elsewhere 
it has been found that it is possible to observe dynamical tun- 
nelling in regimes with much larger atom numbers l35ll . 



V. CONCLUSIONS 

We have studied dynamical tunnelling in a driven ID 
single-well potential provided by the magnetic field of a 
current-carrying wire on an atom chip. An experiment using 
such a geometry would provide a stronger case for dynami- 
cal tunnelling than those performed in optical lattices 10, [H], 
where some ambiguity exists as to the single-well nature of 
the dynamical tunnelling and the correspondence to the clas- 
sical dynamics. 

We have analyzed the dependence of the dynamical tun- 
nelling rate with the experimental parameters and observed 
that it fluctuates considerably, a feature previously observed 
for other trapping potentials. In the limiting case of the quan- 
tum/regular regime, we were able to relate the variation in the 
tunnelling rate to the energy spectrum of an integrable Hamil- 
tonian. In the classical/chaotic regime, the variation in the 
tunnelling rate is linked to a chaos-assisted tunnelling phe- 
nomenon for high H e ff. For lower S e ff we have underlined the 
role of classical resonances linking the islands of regular mo- 
tion with the sea of chaos. 

Finally we have considered the experimental prospects for 
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realising dynamical tunnelling with a BEC in a magnetic mi- 
crotrap. We have found that dynamical tunnelling would be 
observable on a 10 ms timescale for realistic atom-chip pa- 
rameters provided the dynamics in the plane perpendicular to 
the motion could be frozen out. However, this requires BECs 
with a rather small number of atoms. Our conclusions hold 
in the presence of nonlinearities due to atomic interactions as 
long as these are not too large. 
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Appendix A: Atom-chip Hamiltonians 



where g is the Lande g-factor, m F is the magnetic sublevel that 
is trapped and hb is the Bohr magneton. For (B'/Bo) 2 r 2 1 
a Taylor expansion gives 



B 

V(r) = -gm F ^ B —r 2 + const.. 
I Bq 



(A6) 



i.e. the minimum of the potential is approximately harmonic 
with a radial trapping frequency of 



, l gm F ji B 



(A7) 



CO r = B' 

The system Hamiltonian is then 

2 

H = i- + gm Ff x B B V(B'/B ) 2 r 2 + 1. (A8) 
2m 

We define the dimensionless quantities 



B' 



B' - B' 2 

P> H = — —H 



f — — r, t = Q.t, p 

Bo BomQ, mii^ifg 

(A9) 

where £2 will be the angular frequency of the modulation. Af- 
ter dropping tildes, we find our dimensionless Hamiltonian is 



We consider atoms trapped in a magnetic potential formed 
by a current flowing in the positive z direction along single in- 
finite conductor, with a constant bias field B\, in the negative y 
direction. This forms a 2D potential in the x and y dimensions. 
For now we will ignore any dynamics in the z dimension — 
this is reasonable as on typical atom chips the trapping poten- 
tial in this dimension is much weaker and the dynamics are 
correspondingly slower. 

The magnetic field strength a distance yjx 2 + y 2 from the 
conductor will be 



\B(r)\ = 



Mo 



2n 



(Al) 



The bias field B/, cancels the field from the wire along the line 
defined by 



x - x = 



Vol 
2irB h ' 



where the field gradient will be 



\B'\ 



2nxl 



Bb 
xq ' 



(A2) 



(A3) 



We define x' = x - xq and r = a/x' 2 + y 2 . If we add a 
small, possibly time-dependent offset field in the z dimension 
of magnitude Bo then the magnitude of the field at a distance 
r from the minimum can be approximated by 



B 



2 + ^ 



(A4) 



which will be valid for r <K xq (for more detail see the next 
section). The radial potential that the atoms see will be 



V(r) = gm F fi B B V(B'/B )V + 1, 



(A5) 



H=?-+k Vr 2 + 1, 
2 

where we have defined the dimensionless parameter 

_ gm F yi B B' 2 _ ar^ 
K ~ mB Q. 2 ~ Q 2 ' 

The effective Planck's constant for this system is 

hB' 2 n hco r 
«eff = ~ - = yx- 



BlmQ. gm F p B B ' 



(A10) 



(All) 



(A12) 



1. Modulation of the trapping potential 

We now consider the possible ways we can modulate the 
trapping potential for the atoms. The magnetic field at loca- 
tion (x, y) is given by 

B^-^M) - Bb(t) , Bo (tl (A13) 

yjx 2 + y 2 yjx 2 + y 2 J 

where Bo(f) is the offset field along the z-direction, B w {t) is the 
field from the ID conductor, and B^(f) is the bias field along 
the y-direction. We have 



B fe 0D = B fc0 (l +f h \ 
B (f) = B (l+/ Z ), 



2n yjx 2 + y 2 
XoBbo 



^x 2 + y : 



(1+fw), 



(A14) 
(A15) 



(A16) 
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where we have used the shorthand f a = f a (f) for any modula- 
tion for t > 0. We also define the constant 



2nB M ' 



(A17) 



which is the location of the trap minimum with no driving. 
The magnitude of the magnetic field can be shown to be 



B 2 = Bl(\+f z y + 



B lo 



(xq + x') 2 + y 2 



[(i + /*)V 2 + /) 



+ (fb - f w ) 2 4 + 2(1 + f h )(f b - f w )x'x ]. (A18) 
where x' = x - xq. 

2. Choice of time-dependent modulations 

Various Hamiltonians can be realised with appropriate 
choices of modulation of the magnetic fields and currents. 
Here we detail two that are potentially of interest. 

1. fb=f w = fz- 

In this case the trap minimum is stationary, and both 
the field gradient and the offset field Bo are modulated. 
This yields the magnetic field 



B = B (l+fb)(l+r 2 ) l/2 , 



(A19) 



which has a similar form to the modulated standing 
wave used in Refs. J3,Htl — a potential with a stationary 
minimum but a modulated strength. 

2. f b = 0,f z =0: 

This case is probably the easiest experimentally, as only 
the current in the trapping wire needs to be modulated. 
This avoids any difficulties with inductances in the coils 
providing the bias and offset magnetic fields. The field 
in this case is 



B = B °{ 1+ (S [y 2 + ( x '- x °f^ 2 } 



1/2 



(A20) 



This potential has a constant strength but the position of 
its minimum oscillates. 

In summary, the experimentally relevant Hamiltonians in 
dimensionless units are 

H = ^+*(l+ecosT)[l+(x' 2 +3; 2 )] 1/2 , (A21) 



H = — + k [l + y 2 + {x - xoecos r) 2 J 



211/2 



(A22) 



In this paper we have only considered Hamiltonian (IA21I) as 
this is analogous to the experiments performed in the modu- 
lated standing wave. 

Appendix B: Dimensional reduction of the Gross-Pitaevskii 
equation 

We begin with the standard Gross-Pitaevskii equation for 
the mean-field of a condensate tightly confined along the y- 



dimension [3 

d 

+ T2d|(A(jc, z)| 2 ]*A(x, z) 



ti 1 Id 2 d 2 , 



(Bl) 



V(x,z, t) = V [l + ecos(Qf)] Vl + (x/d) 2 , (B2) 



with 



The effective 2D interaction is j2D = 2 ~\j2nti 2 a s / '(ma y ) Q4; 
where a s is the 3D scattering length and a y = ^JWJJnujjyj the 
oscillator length in the frozen y-direction. J dxdz\4>\ 2 = N, the 
total 3D number of atoms. The parameters of the potential can 
be read from Eq. dA8b with the identification x <-> r. 

We now measure length in units of d, time in units of Q. 1 
[see Eq. ([3])] and employ the dimensionless wave function tfr = 
kjj y/d/N, which is normalized to one. Finally relabelling tfr — > 
kjj we obtain: 



d 



[4 


f- 


d 2 \ 




[dx 2 


+ dT 2 ) 



where 



+ U 2D \Hx,z)\ 2 ]t/f, 
V{x,z,i) = k[\ + ecos(0] Vl + x 2 . 



(B3) 



(B4) 



Importantly U2D = Ny2D/mQ 2 d 4 . We now eliminate the di- 
mension z, assuming the size of the condensate in that direc- 
tion, L z , to be large and the dynamics hence slow. We simply 

set \f/(x, z) = tff(x)4>(z), with <f>{z) = 6(L/2 - \z\)/ VL, where 6 is 
the Heaviside function and L = L~/d. Insertion into Eq. (IB3b 
yields: 



+ d - 
ot 



d 2 



-^f^ + V(x,t) + U lD \{(r(x)\ 2 
2 ox 1 



h (B5) 



with [/id = U2d/L z . We obtain Eq. (fT4t by replacing x — > q 
and have 



U 



ID 



Nj 2 D 2 ^h 2 a s N 2 yfhiti 2 a s NK 



mQ. 2 d 3 L z m 2 a v Q. 2 d 3 L z ma r dL z V u 



(B6) 



To rewrite Eq. ( IB6b in terms of the quantities h e g, k, e and 
U\d that underly most of our numerical results, we substitute 
d = Bq/B', Vq = gmffisBo [comparing Eq. (IB2b andEq. ( IA8H 
and then use Eq. dA12b to arrive at: 



L z U m 



2^hia s h 1 'J 2 K x ^ \ 



(B7) 



The atom numbers and densities quoted in section HV Cl corre- 
spond to L z = 50 fim, to, — 2n x 18.6 Hz, u y — 2n x 30 kHz 
and other parameters as given in section HV Bl 
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